Impurity and correlation effects on transport in one-dimensional quantum wires 
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We study transport through a one-dimensional quantum wire of correlated fermions connected to 
semi-infinite leads. The wire contains either a single impurity or two barriers, the latter allowing 
for resonant tunneling. In the leads the fermions are assumed to be non-interacting. The wire is 
described by a microscopic lattice model. Using the functional renormalization group we calculate 
the linear conductance for wires of mesoscopic length and for all relevant temperature scales. For a 
single impurity, either strong or weak, we find power-law behavior as a function of temperature. In 
addition, we can describe the complete crossover from the weak- to the strong-impurity limit. For 
two barriers, depending on the parameters of the enclosed quantum dot, we find temperature regimes 
in which the conductance follows power-laws with "universal" exponents as well as non-universal 
behavior. Our approach leads to a comprehensive picture of resonant tunneling. We compare our 
results with those of alternative approaches. 



I. INTRODUCTION 

The interplay of static impurities and correlations in 
one-dimensional (Id) Fermi systems leads to a variety 
of surprising effects. A detailed understanding of the 
physics involved is interesting from the basic-science 
point of view as well as from the perspective of a pos- 
sible application of Id quantum wires in future nano- 
electronics. In a Id metallic system, correlations have a 
strong effect on the low-energy properties. Quite differ- 
ently from the conventional Fermi-liquid behavior and 
even in the homogeneous case, physical properties at 
low energy scales follow power-laws, described by the 
Tomonaga-Luttinger liquid (TLL) phenomenology.^ For 
spin-rotational invariant and spinless models the expo- 
nents can be expressed in terms of a single number, the 
interaction-dependent TLL parameter K. 

A single static impurity in a TLL with repulsive in- 
teraction (Q < K < 1) leads to a dramatic modification 
of the low-energy physics as is obvious from lowest order 
perturbation theory in the strength of the impurity. 
More elaborate methods show that for vanishing energy 
scale s the conductance of such a system vanishes follow- 
ing a power-law. The scaling exponent 2aB is inde- 
pendent of the bare impurity strength and determined 
by the exponent as = 1/K — 1 (for spinless fermions) 
of the one-particle spectral weight of a TLL close to an 
open boundary.^ For weak impurities the correction to 
the impurity free conductance scales as s'^^-'^~^\ which 
holds as long as the correction stays small. Assuming 
an infinite system the asymptotic low-energy properties 
have been investigated intensively within an integrable 
field theoretical model — the local sine-Gordon model 
(LSGM).^"^'i° For this model, with fixed K the conduc- 
tance as a function of temperature follows one-parameter 
scaling for different strengths of the impurity. Similar re- 
sults were obtained for a TLL connected by arbitrarily 
"smooth" (spatially adiabatic) contacts to semi-infinite 
Fermi-liquid leads. ^^"^^ 

Also transport through a double barrier enclosing a 



quantum dot has been studied over the last few years 
using field theoretical models. ^'^^"^^ Tuning the dot en- 
ergies by a gate voltage Vg resonant tunneling can 
be achieved. For appropriately chosen dot parameters 
(dot size, barrier height) and as a function of temper- 
ature T the conductance G shows different temperature 
regimes with distinctive power-law scaling and exponents 
which can be expressed in terms of ^.17,19-21,24-27 p^j. 
the double-barrier problem no exactly solvable generic 
model is known and applying different approximate 
analytical^^d9-2i,24,25,27 g^^^^ numerical^^ methods has 

not provided a consistent picture. 

To study the effect of the interplay of correlations and 
impurities on transport on all relevant energy scales we 
use the functional renormalization group (fRG) method. 
The fRG can be applied directly to microscopic models. 
We consider the lattice model of spinless fermions with 
nearest-neighbor hopping t and nearest-neighbor inter- 
action U. The impurities are modeled by either locally 
raising site energies or reducing hopping matrix elements 
across bonds. The fRG was recently introduced^^ as a 
new powerful tool for studying interacting Fermi systems. 
It provides a systematic way of rcsumming competing 
instabilities^^ and goes beyond simple perturbation the- 
ory even in problems which are not plagued by infrared 
divergences. The fRG procedure we use starts from an 
exact hierarchy of differential flow equations for the one- 
particle irreducible vertex functions, ^^"^'^ as e.g. the self- 
energy and the effective two-particle interaction. It is 
derived by replacing the free propagator by a propagator 
depending on an infrared cutoff A and taking the deriva- 
tive of the generating functional with respect to A. 

At T = we introduced earlier two truncation 
schemes which led to a manageable number of coupled 



equations 



34-36 



The flow of the self-energy, which in par- 



ticular encodes the renormalized impurity potential, is 
fully taken into account, while the two-particle vertex is 
parametrized by a single flowing coupling constant. The 
bare two-particle interaction U is taken as a small param- 
eter, but the impurities can have arbitrary strength and 
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shape. We are thus in a position to study the single- and 
double-barrier problem by applying the same approxi- 
mate method. For weak and strong single impurities the 
flow equations can be solved analytically."^^ The results 
are consistent with the above mentioned power-law scal- 
ing of the conductance. For general parameters, in par- 
ticular in the double barrier case, we solved the coupled 
differential equations numerically which at T = can eas- 
ily be done for systems of up to = lO'' lattice sites. Us- 
ing the fRG we previously studied the one-particle spec- 
tral function close to a single impurity,^*"^^ the decay 
of density (Friedel-) oscillations off an impurity,"^^ and 
the persistent current in mesoscopic rings pierced by a 
magnetic flux.^'^ For these observables we recovered the 
expected "universal" power-laws with exponents which 
for 1/2 < K < 1 {K ^ 1 ioT U ^ 0) turned out to be in 
good agreement with the ones known from the LSGM. In 
addition it was shown that the asymptotic behavior typ- 
ically sets in on very small energy scales and for very 
large systems, except for very strong bare impurities. 
Within our approximation, in the calculation of various 
observables we have to solve as a last step the problem 
of a non-interacting fermion moving in the renormalized 
impurity potential. This mapping on an effective one- 
particle problem is helpful for the understanding of the 
observed effects. 

In a short publication we previously discussed the 
transport for the single-barrier case.^^ The results were 
obtained after supplementing the A'^-site wire of interact- 
ing fermions (called the "interacting wire" in the follow- 
ing) by Id semi-infinite leads. We studied the T = 
behavior as a function of For "smooth" contacts 

between the leads and the interacting wire our data for 
different impurity strengths and inverse system sizes col- 
lapse onto a single curve after applying a one-parameter 
scaling ansatz. This shows that even a small bare im- 
purity leads to a vanishing conductance in the limit 
N —i oo. Kt K = 1/2 the one-parameter scaling func- 
tion of the infinite LSGM is known analytically.^' ^'^ Up 
to a small error our rescaled data fall onto this curve. 
As an additional benchmark for small wc compared 
the fRG approximation for the conductance with re- 
sults of density-matrix renormalization group (DMRG) 
calculations. ^^'^^ These tests lead us to conclude that our 
approximate method works reliably for interactions in the 
range l/2<if<l. In addition, the collapse of our fRG 
data onto the K = 1/2 local sine-Gordon scaling function 
provides an indication that the low-energy physics of the 
LSGM can indeed be found in a large class of models of 
inhomogeneous, one-dimensional, and correlated electron 
systems, which is not obvious a priori. 

In the present publication we extend the fRG scheme 
of Ref. 36 to finite temperatures. As in any experimen- 
tal system the bandwidth B and the inverse of the wire 
length set upper and lower bounds for the possible "uni- 
versal" scaling behavior of the conductance. This has to 
be contrasted to studies of field-theoretical models with 
B = oo and (in most cases) 1/N = 0. We calculate G 



on all temperature scales. Depending on the bare impu- 
rity and interaction strength an asymptotic low-energy 
regime might not be reachable in experiments on finite 
wires. For finite T we consider up to A?^ = 10^ sites. For 
typical lattice constants this leads to interacting wires 
in the micrometer range, roughly corresponding to quan- 
tum wires accessible to transport experiments. Within 
the last few years, transport properties were measured in 
different, effectively Id mesoscopic wires. 

Although in many respects our model is more realistic 
than the ficld-thcorctical models studied so far, impor- 
tant ingredients for an appropriate description of the ex- 
perimental situation are missing, as e.g. the spin degree 
of freedom, higher-dimensional leads, and realistic con- 
tacts. In the present paper we therefore refrain from a 
detailed comparison of our results to experiments. Later, 
we briefly comment on certain aspects of the experimen- 
tal results and discuss the possibility of extending our 
method to treat models closer to experimental systems. 

In the two limits of weak and strong impurities we re- 
cover the expected power-law scaling of G{T) with expo- 
nents which are in good agreement with the ones known 
from the LSGM. We show that one-parameter scaling is 
achieved. We also present results for attractive interac- 
tions with U < and K > 1. Based on our findings 
we discuss the accuracy of our approximation scheme in 
detail. We investigate the relation between our approach 
and a "leading-log" resummation for the effective trans- 
mission at the chemical potential introduced in Ref. 9. 
We then consider the double-barrier problem. Within 
a single approximation scheme all relevant temperature 
regimes can be investigated, leading to a comprehensive 
picture. We compare our findings to results that have 
been obtained for restricted temperature regimes by ap- 
plying alternative approximate techniques and a numeri- 
cal method to field- theoretical models. A brief account of 
results for resonant tunneling in the double-barrier prob- 
lem obtained from the fRG has been presented earlier, 
focusing on the scaling of the peak conductance for sym- 
metric barriers. 

This paper is organized as follows. In Sect. II we intro- 
duce the model considered. In Sect. Ill results of single- 
particle scattering theory are derived. They are useful for 
the understanding of our findings for the conductance. 
The details of the fRG method are presented in Sect. 
IV. We discuss the general approximation scheme and 
describe how the fRG is extended to finite T and to the 
study of transport. In Sect. V we present our results for 
the transport in the presence of a single impurity as well 
as through a double barrier. We conclude with a sum- 
mary and an outlook in Sect. VI. In the appendix we 
briefly describe the relation between the imaginary part 
of the self-energy and the occurrence of current-vertex 
corrections in a Landauer-Biittiker-type expression for 
the conductance. 
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II. THE MODEL 

We consider spinless fermions on a lattice with nearest- 
neighbor interaction supplemented by various types of 
impurities. The Hamiltonian reads 

H = ifkin + -f^int + -f^imp • (1) 

The kinetic energy is given by nearest-neighbor hopping 
with an amplitude i, 

oo 
j=-oo 

where we used standard second-quantized notation with 
and Cj- being creation and annihilation operators on 

site j, respectively. 

The part of the Hamiltonian containing the interaction 

reads 

Hint = J2 ["'J ~ ~ 

with the local density operator rij = Cj Cj . The inter- 
action acts only between the bonds of the sites 1 to N, 
which define the interacting wire. The interaction Ujj+i 
between electrons on sites j and j + l is allowed to depend 
on the position. As discussed in Refs. 38 and 39 in the 
absence of a one-particle impurity potential, turning on 
the interaction sharply leads to a reduction of the con- 
ductance compared to e'^/h.*^ As in the other theoretical 
studies mentioned we here avoid this effect by modeling 
"perfect" contacts. To achieve this Ujj+i is taken as a 
smoothly increasing function of j starting form zero at 
the bond (1,2) and approaching a constant bulk value U 
over a sufficiently large number of bonds. Equally, the 
Ujj+i are switched off close to the bond {N — 1, N). The 
results are independent of the detailed shape of the en- 
velope function as long as it is sufficiently smooth. The 
larger N the "smoother" Ujj+i has to be chosen. More 
explicitly we here take {N even, j = 1,. . . , N/2) 

jj ^ jj arctan [s{j - js)] - arctan [5(1 - j^)] 
■''■'"^^ arctan [s{N/2 - j^)] - arctan - js)] 

for the left part of the interacting wire and likewise for 
the right. For odd N a similar function is used. The 
parameters s and js are chosen such that in the absence 
of any impurity for T = we find 1 - G/{e^/h) < 10"^. 
To achieve this for the A'' considered here js = 56 and 
s « 1/4 is sufficient. The two regions of the lattice with 
j < 1 and j > N constitute the two semi-infinite leads. 

In Eq. (3) we shifted rij by a parameter i'{n, U), which 
depends on the filling factor n and U. This is equivalent 
to introducing an additional one-particle potential in the 
interacting part of the system. The parameter u{n, U) is 
chosen in such a way that the density in the interacting 
wire acquires the desired value n. 




FIG. 1: Schematic plot of the quantum dot situation, where 
the barriers are modeled by two site impurities. 

The general form of the impurity part of the Hamilto- 
nian is written as 

hi' 

where Vjj' is a static potential. Site impurities are given 
by a local potential 

Vj.r = V, 5,,,. . (6) 

In our study of the effect of a single barrier we mainly 
consider a site impurity located far away from both leads 
on site jo with jo » 1 and N — jo ;:g> 1, 

Vj = V5jj,. (7) 

In the discussion of resonant tunneling two site impurities 
of strengths Vi and Vr on the sites ji = ji — I and j,. = 
jr + 1 constitute the barriers, with ji 1 and N—jr^ 1. 
The Nd sites between ji and jr define the quantum dot 
with jr — ji + 1 = Nz) > 1. The effect of a gate voltage is 
described by a constant Vg on sites j; to jr- This situation 
is depicted in Fig. 1. We also consider hopping impurities 
described by non-local potentials 

as barriers. For the special case of a single hopping im- 
purity, 

= - *) ^jJo > (9) 

the hopping amplitude t is replaced by t' on the bond 

linking the sites jo and jo + 1. In the double-barrier prob- 
lem we consider a hopping across the bond {ji , ji ) and 
tr across (jnjr)- For small ti and tr also the interactions 
across the bonds {ji,ji) and (jr, jr), i-e. Uj^ j^ and Uj^j^, 
are assumed to be reduced. In the following we set the 
bulk hopping amplitude t equal to one, i.e. all energies 
are given in units of t. This leads to the non-interacting 
dispersion Sk = — 2cosfc and bandwidth B = 4. The 
lattice constant is set to 1. 

The homogeneous model H = if kin + -ffint with a con- 
stant interaction U across all bonds (not only the ones 
within [1,A'']) can be solved exactly by Bethe ansatz.^° 
It shows TLL behavior for all particle densities n and 



4 



any interaction strength except at half fiUing for \U\ > 2. 
The U - and n-dependent TLL parameter K can be deter- 
mined solving coupled integral equations, ''^ which in the 
half-filled case can be done analytically with the result 



K = 



arccos 



(10) 



for |J7| < 2. 

We compute the conductance of the above models in 
the linear-response regime. For T = the fRG flow equa- 
tions can be solved numerically for interacting wires as 
large as = 10'' lattice sites. For finite T we treat sys- 
tems with up to A'' = 10* sites. 



III. SCATTERING ON ONE-DIMENSIONAL 
LATTICES 

The truncation of our fRG flow equations leads to a 
scattering problem of non-interacting fermions in an ef- 
fective potential given by the non-trivial spatial depen- 
dence of the self-energy on sites 1 to jV. 34-36 j|- jg t^^us 
useful to present some results for the scattering in a Id 
system with an arbitrary (non-local) potential in a finite 
segment of the infinite wire. The results are important for 
the understanding of how the fRG based approximation 
scheme is set up and how the conductance is calculated 
within this scheme (see Sect. IV). In the first subsection 
we focus on this aspect. In the second more technical 
subsection we derive relations which can be used to gain 
a deeper understanding of our results for the conduc- 
tance presented in Sect. V. In order to make this section 
self-contained we include short derivations of well-known 
results, as e.g. Eq. (15). 



A. General relations 

The single-particle Hamiltonian Hip = Hq + Vlr we 
consider in this section reads 



= - E (li-i)01 + ii.c.) 



3=-oo 
oo 



- E (li)0- + i| + ii.c.) + if. 

Flr = -fL(|0)(l|+h.C.)-ffl(|iV)(iV + l|+h.C.) , 

(11) 

where \j) denotes the state in which the fermion is cen- 
tered at lattice site j. Here Hs is a general single-particle 
Hamiltonian in the scattering segment from site 1 to 
which is connected to ideal leads with nearest-neighbor 
hopping t = 1. The hopping matrix elements II and Ir 
connecting the three parts described by Hq are assumed 
to be real. 



There are two ways to describe the scattering problem. 
In the usual approach the ideal infinite system is used as 
the unperturbed system and the deviations are treated 
as the perturbation. In Eq. (11) the grouping has been 
made differently. The disconnected system "left lead - 
scattering region - right lead" is described by and 
the connections proportional to II and Ir present the 
perturbation in the construction of the scattering states 
|A;,a±)5i 



|fc,a±) = lim 



|A;, a) 



v^o Ek - Hip ± irj 
?(efc ± IT]) 

k,a) +G{sk±iO)VLR\k,a} . 



= lim ztirj G{ekii'r])\k, a) 



(12) 



Here the index a stands for L, R and the unperturbed 
states \k,L{R)) are standing waves in the left (right) 
semi-infinite lead, e.g. for j <0 



(13) 



{j\k,L} = d-sm[k{j-l)] 



In this section G{z) denotes the resolvent matrix with 
respect to Hip. With the above normalization of the 

standing waves, {k,a\k',a') = 5a.a'5{k — k') holds. For 
j > N the scattering state \k, L+) has components 



{j\k,L+) = ^{j\G{ek + tO)\l)tL{0\k,L) 
= {j\Goiek+iO)\N + l)tii 
x{N\G{ek+iO)\l)tL{0\k,L). 



(14) 



In the second equation we used G = Go + Go Vlr G with 
Go{z) = {z — Ho)~^. The resolvent matrix element of the 
semi-infinite chain in Eq. (14) is given by^^ 



{j\Go{ek+iO)\N + l) 



ikU-N) 



which shows that {j\k,L+) is in fact an outgoing plane 
wave. From Eqs. (13) and (14) the transmission ampli- 
tude t{£k) can be read off. The corresponding probability 
is given by 



|i(efc)|' =4ii4sin^fc |(7V|G(efe + z0)|l)r 



(15) 



This is a very important relation, as in the Landauer- 
Biittiker formula^'^ the transmission probability di- 
rectly determines the linear conductance G(T) of non- 
interacting fermions at temperature T. For spinless 
fermions on the lattice it reads 



G(T) 



B/2 
B/2 



de 



\t{e)\^de, 



(16) 



where /(e) = l/{e^^^~'^^ + 1) is the Fermi function with 

P = 1/T and ^ the chemical potential. The Boltzmann 
constant ks is set to 1. If the fRG fiow is integrated for 
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a finite temperature and a finite interacting wire the self- 
energy is T- and A^-depcndent, i.e. and therefore also 
the transmission probability depends on tempera- 

ture and wire length. In the notation used in this section 
these additional dependences are left implicit. 

For the further discussion of the resolvent matrix ele- 
ment in Eq. (15) the relation^^ 



PG{z)P = 
zP 



(17) 



PH^pP ■ 



PH^pQ {zQ - QHipQ)-^ QH^pP 



where P and Q are projectors with P + Q = 1, is used 
several times. With this projection formula the calcula- 
tion of {N\G{ek + iO)|l) for the infinite system can be 
reduced to the problem of inverting an TV x A''-matrix. 



With the definition P, = J2j=i li)01 

one obtains 

{N\G{z)\l) = {N\ [zPs - Hf{z)]-' |1), (18) 

where is the diagonal element of the resolvent of the 
semi-infinite chain at the boundary and 

Hf{z) = H, + Gtiz) {tlm\+tl\N){N\) . 
For energies inside the band G° is given by 



Gg(£ + iO) = -e'^^'^ = {£- iV4-£2)/2 . 

In our fRG description of interacting fermions the 

"scattering segment" [1, N] is defined by the part of the 
wire where interactions are present. 



B. The double barrier problem 

As one of the problems in our paper we discuss trans- 
port through a quantum dot with No sites extending 
from j; ^ 1 to with N—j,. 1, where jV = Ji+Nd — I. 
For the sites neighboring the dot we use the indices ja 
(o = l,r). In Hs the (real) hopping matrix elements 
from site ji to ); are denoted by t; and from jV to jV by 
fj.-^^ Proceeding similarly as in Eq. (14) one can express 
Gat,! = (iV|G|l) as 



Gn,i 
where 



-GN,j,iiG^.'\ 



G^;.>G,,,,,t;G°;;^, (19) 



G°'''{z) = [z-H^,piia=0)]- 



be reduced to the inversion of an No x Nd matrix. With 
Pdot = \j)U\ one obtains 

PdotG(^)Pdot = [zPdot - PdotHsPdot - Hbi^T' , (20) 



where the boundary term is given by 

Ht{z) = Ti{z)\ji){ji\+T,{z)\jr){jr\, (21) 

with 

r<,(z) = M';, (z). 

Ja^Ja 

As the information from outside the dot is contained 
in Fa it is desirable to express the conductance also in 
terms of these functions. If one defines the projector 
Pi = Sj=i using Eq. (18) one can derive the 

following identity for G^^p^ = PiG°^^Pi, 

G?^V,(^ + ^0)-G?^P,(^-^0) = 
-ziiv/4^G°^ {s + iO)|l)(l|G?;,'p, {e - i0).(22) 



The matrix element of this equation relates |G' 



0,1 |2 



to the imaginary part of G^^' . With the analogous rela- 
tion for the right neighbor of the dot one finally obtains 
using Eqs. (15) and (19) 

\t{s)f=4At{s)Ar{e)\Gj^,j,{s + iO)f. (23) 

Here we have decomposed 

with fio and Aa real functions. This equation is a useful 
generalization of Eq. (15). 

The explicit result for the transmission probability is 
simplest for the case Nd = 1. Then ifdot = -Pdot -f^s -Pdot 
consists of a single term^^ and one obtains the 

generalized Breit-Wigner form^^ 



4Ai{s)A,{e) 



e-Vg-ni{s)-nr{e] 



is the resolvent for the system cut into two pieces by 
putting ta = 0. Since G'^'" enters the expressions for 
the transmission probability derived below, we emphasize 
that G°'° is not the Green function for the Hamiltonian 
(11) but an auxiliary Green function. Using the projec- 
tion technique again the problem of calculating G^^j, can ^ _ 



+ [Ai{s) + Ar.{e)f 

(24) 

In all cases discussed in this paper a resonance occurs at 
the energy en, where en — Vg — 0,i{eii) — flri^R) vanishes. 
Only for a symmetric dot with A/ = = A the peak 
value is given by one, i.e. perfect transmission. Later 
Eq. (24) will also be used for the understanding of the 
conductance in the presence of a single site impurity. 

For general dot size, Gj^j^ can be expressed in terms 
of the Ta and the resolvent of the isolated dot, 

G°iz) = {zPiot - Haot)-^ , 
using the special form of Hi, in Eq. (21) 



(1 - r^GO ,,,)(i - r,Go^,,j - r.r.GO ,,^G^,,, ' 

(25) 

For |r;|,|rr| ^ 1 one can have Njj narrow resonances 
which have the generalized Breit-Wigner form as in 
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Eq. (24). This follows from Eq. (25) by using the spectral 
representation 



GlAz) 



(j|g«)(ga|/) 



(26) 



where the \ea) arc the cigenstatcs of iJdot- The corre- 
sponding eigenvalues depend linearly on Vg. For dots 
weakly coupled to the rest of the system the poles in 
Eq. (26) lead to narrow resonances in Gj^.,j,. For ener- 
gies close to Ea the resolvent matrix elements G'jj, can 
be replaced by the corresponding single-pole term. Then 
the terms quadratic in (5° in the denominator of Eq. (25) 
cancel and for e one obtains 



(27) 



4A|")(£)Ar^(e) 



Ml 



£-£„-^^"^(£)-^r(£) 



A["^(£) 



where A^\s) = |(ja|£a)pAa(£) and correspondingly for 

rii"''. The factor |(.ja|£a)P is typically of order 1/No 
and the widths of the resonances are reduced by this 
factor compared to the No = 1 dot. As the separation 
of the resonances is of order B/Njj the generalized Breit- 
Wigner form as in Eq. (24) holds for [Fil, IF^I < 1. This 
can be fulfilled for <^ 1 or for large site impurities 
at sites ji and jr (see below). If ra(£) is assumed to be 
slowly varying over energies of the resonance width the 
integrated weight Wa of the Lorentzian resonance near 
Ea is given by 



t«„«47r[l/Aj"'(£,) + l/A(«)(£„)]- 



Here we have assumed that 1 — (fio 



1, which holds 



for large dots with separated resonances. 

Using the general projection formula (17) a relation 
between the diagonal matrix elements of the auxiliary 
Green function G?'°^ to the left and right of the dot 

[entering Eqs. (24) and (27)] and the corresponding ma- 
trix elements of the exact Green function 5 can be 

_ JatJa 

derived. For general dot size and with a one obtains 



(1 - T/GO ,J (1 - TrGlJ - r,r.GO ,,^G5^,,, 





£-£„-Oi"H£)+iAi")(£)' 


G?'". (£ + iO) 


e - Ea - ^'i"\s) - nl"\e) + I 


A["'(£)+A[")(£)_ 



(28) 

where the second line holds tov z = e + iO Sa in the 
case of dots which are only weakly coupled to the rest of 
the system. 

For the interpretation of the temperature dependence 

of the conductance of a dot with narrow resonances in 
the regime B/Nd <^ T <^ B also the first equality in 



Eq. (19) is useful. Applying the general projection for- 
mula (17) yields 



G 



0,1 

N,]i 



1 - ifG°:'.. G°f. 

' 31, Jl 



With Eq. (22) and the corresponding relation for G^'^^ 
one obtains another exact relation for the transmission 
probability 



4A,(£)ImG°;^(£ + iO) 



l-ri{e + iO)G°fj^{e + iO) 



(29) 



If IF; I <C 1 but ir^l w 1 no resonances occur and pertur- 
bation theory in Fj is possible. Putting the denominator 
in Eq. (29) equal to 1 leads to the "golden-rule" approxi- 
mation. The conductance for this "single-impurity case" 
is then to leading order in if given by 



(30) 



f-S/2 
3/2 

xImG°;^ (£ + ^0)d£. 



We now return to the dot situation with narrow reso- 
nances. Its simplest realization is given by two hopping 
impurities with tf,t^<^lm an otherwise perfect lattice. 
Then the G?'% and \{ja\£a)\^ ioi a = l,r are identi- 
cal and the expression for the weights of the resonances 
simplifies to 



Wo 



Ki;|£„)|2|ImG°;' (£„ 



For B/Nd < T <C B the transmission probability |^(£)P 
in the Landauer-Biittiker formula (16) can approximately 
be replaced by J2a "^a^i^ ~ This yields 



5dot(T) 



1 1 

f2 + ^ 



B/2 
B/2 



xImGlj^{e + iO)de. 



-f)lraG':'. (s + zO) 
de J 3i,jr 

(31) 



The integrands in Eqs. (30) and (31) differ in the last fac- 
tor. For the simple example with ^ 1 in an otherwise 
perfect lattice G^''^-^ is just the boundary Green function 
G° of a semi-infinite lattice. For Nd 3> 1 and i ^ <C 1 the 
imaginary part of the boundary Green function Im G"^ 
of the dot can for the purpose of integration in Eq. (31) 
be approximately replaced by ImG°. In Sect. V we en- 
counter more interesting realizations for which the sim- 
plified expression for Wa as well as this replacement can 
be used. Comparison of Eqs. (30) and (31) then shows 
that in the "high-temperature regime" B/Nd -^T B 
the conductance of the dot can be obtained by adding 
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resistances for the single-impurity problems (Kirchhoff 's 
law) described by Eq. (30), 



(32) 



<7dot(r) 9i{T) gAT)' 



The quantum mechanical interference effects in the dot 
are washed out for large enough temperatures. For lower 
temperatures the approximation Eq. (32) breaks down. 
Another temperature window where the result for the 
conductance takes a simple form is given by | Ai"^ (£„)] ^ 
T <C B/Nd- If the resonance is at the chemical poten- 
tial or much closer to it than T the conductance decays 
with temperature oc Wa/T. For the interacting case dis- 
cussed in Sect. V the weight of the resonance Wa itself is 
temperature dependent which results in G{T) oc T"'^~^. 

The above arguments for weak hopping connections of 
the dot with the rest of the system can easily be gen- 
eralized to the case of « 1 but a large site impurity 
Va at the site ja- For |V;| 1 the site ji can be "inte- 
grated out" using Eq. (17). This leads to a weak hopping 
^ji-2,ji = *7i-2,ji-i ii/^i between the sites ji-2 and ji . 

The relations presented in this section play an impor- 
tant role for the interpretation of our fRG results for the 
transport through interacting quantum wires. 



IV. THE FUNCTIONAL RENORMALIZATION 
GROUP AT T > 

We now return to our Hamiltonian Eq. (1). Due to 
the presence of leads the direct calculation of the non- 
interacting propagator related to ifkin + ^^imp [see Eqs. 
(2) and (5)] requires the inversion of an infinite matrix. 
Using Eq. (17) it can be reduced to the inversion of an 
N X N matrix. The leads then provide an additional di- 
agonal one-particle potential on sites 1 and N, depending 
on the Matsubara frequency w„, 




(33) 



Since the interaction is only non- vanishing on the bonds 
between the sites 1 to A'' the problem including the semi- 
infinite leads is this way reduced to the problem of an 
A'^-site chain. The chemical potential of the entire system 
in thermal equilibrium is denoted by /i(n, U, T). 

To set up the fRG the projected non-interacting prop- 
agator Go is replaced by 



(34) 



with a function which is unity for \uJn\ 3> A and van- 
ishes for \u>n\ -C A. At T = we earlier considered 
X^{oj) = e(|a;| - A).3^-36 For technical reasons at T > a 
smoother cutoff function, specified below, is more appro- 
priate. We here choose to include the impurity part of the 



Hamiltonian iJimp in G^. Earlier it was taken as a part 
of the self-energy.^® At T = both choices lead to the 
same results. By means of Gq the generating functional 
of the one-particle irreducible vertex functions for the 
Hamiltonian H becomes A dependent. Differentiating 
the functional with respect to A and expanding it in the 
external sources then leads to an exact infinite hierarchy 
of coupled flow equations for the vertex functions. 
In practical applications this set of equations has to be 
truncated. As a first step we neglect the three-particle 
vertex which is small as long as the two-particle vertex 

stays small which is the case for not-too-largc U.^^ 
This approximation leads to a closed set of equations for 

and the self-energy given by 



d_ 
dA 

and 



S'^(l', 1) = - 4 X) ^^^^^ ^' ^) 

^ 2,2' 

(35) 



Ar^(l', 2'; 1,2)= IY.I1 ^^(3' 3') S^{4, 4') 



3,3' 4,4' 



X r^(l',2';3,4)r^H3',4';l,2) 
-r'^(l', 4'; 1, 3) r'^(3', 2'; 4, 2) - (3 ^ 4, 3' ^ 4') 
-hr^(2', 4'; 1, 3) r^(3', 1'; 4, 2) -|- (3 ^ 4, 3' ^ 4') . 

(36) 

The labels 1, 2, etc., stand for the quantum numbers of 
the non-interacting one-particle basis and the Matsubara 
frequencies. The full propagator G^ is given by 



G^= [(G^)-I-EA]"' 
and the so-called single-scale propagator by 



5^ = G^ 



OA 



(37) 



(38) 



1- 



-§aX^(x) 



1/(27iT) 



A-7tT A A+jiT X 



A-jiT A A+jiT X 



FIG. 2: The cutoff function x i^) for fixed A and the negative 
of its derivative (with respect to A) . 

We choose the cutoff function 

0, |a;„|<A-7rT 
X^(a;„) :={ ^ + , A - ttT < K| < A + ttT 

1, \oJn\>A + TrT 
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with the accompanying derivative 

d 



dA 



1 

2-n-T ' 

0, 



A-7rT< 
otherwise , 



< A + ttT 



where A starts at oc' and goes down to 0. Both these func- 
tions are shown in Fig. 2. With this choice and for a fixed 
A a Matsubara sum involving the single-scale propaga- 
tor Eq. (38) contains only the two terms with frequencies 
A — ttT < \ujn\ < A + ttT. This reduces the numerical 
effort in solving the coupled difFercntial equations for the 
self-energy and the two-particle vertex. With this cutoff 
function the right-hand side of the flow equation has a 
jump at every Matsubara frequency and is smooth in be- 
tween, such that the flow can be integrated by standard 
routines. 

To obtain a manageable set of equations even for large 
N we next introduce further approximations in the flow 



equation of the two-particle vertex. We are dealing with 
at most two local impurities. For large N such a small 
number of impurities does not significantly alter the flow 
of the effective interaction and we thus neglect the feed- 
back of the bare impurity potential entering Gq on the 
flow of r^. Self-energy effects on the flow of F^ are of 
order and are ignored. For the integration of Eq. (36) 
we furthermore consider a system with constant U on 
all bonds. The two-particle vertex is projected onto the 
Fermi points and parameterized by a nearest-neighbor 
interaction U^. After taking the limit N ^ oo the flow 
equation for reads^^ 



_d_ 
dA 



A\2 



with 



(39) 



J 



h{x) 



Re 

2n 



- {no + ix) 



'1 - 



(/zo -I- ixY 



3i/x^ - IQnlx - 12ifil{x'^ + 1) + Gx^fJ'O + ISx/Uq -|- 6ix'^ + ix'^ 
7r(2/xo -I- «a:)(4 - nl+x'^ - 2ixiJLQY 



-(40) 



where a)„ stands for the fermionic Matsubara frequency 
closest to A, and /ig = — 2cosA:i?. The above simpli- 
fied fiow equation yields not only the correct low energy 
asymptotics to second order in the renormalized vertex, 
but contains also all non-universal second order correc- 
tions to the vertex at ±kF at higher energy scales. In 
the case where the interaction depends on position, as 
an additional approximation we apply Eq. (39) locally 
for each bond. As long as the bulk part of the interact- 
ing wire is much larger than the contact regions this has 
a small effect. 

Because of the above approximations the self-energy is 
a frequency-independent tridiagonal matrix in real space: 
only Yij j and Y,j j±i are non-zero. In the exact solution 
a frequency dependence of S is generated to order C/^ 
(bulk TLL behavior), which shows that in the approx- 
imation for the self-energy terms of order C/^ are only 
partly included. For the flow of the T- and A'^-dependent 
self- energy we obtain 



d 



St, = - 



dA ^-^ 2tt 



|w„|RiA r=±l 

xQ{iuin) ■ 



Q(iw„) - x^K)SA 



j+r,j+r 



dA i'^^^ 2n 



|a;„|«iA 

xQ{iuJn) 



Q{iuJn) - X^K)SA 

1 



(41) 

where |a;„| w A stands for taking the positive as well as 
negative frequency with absolute value closest to A, and 
J, j±l e [1, A^]- The matrix Q is the inverse of the cutoff- 
independent, projected, and non-interacting propagator 
Q = {Gq)-\ 

The initial condition of the self-energy at A = oo is 
given by 



-v{n, U) {Uj-i,j + Ujj+i) 



= 0, 



and the one of the flowing nearest-neighbor interaction 
by U°° = U. In a numerical solution the flow starts at a 
large finite initial cutoff Aq. One has to take into account 
that, due to the slow decay of the right-hand side of the 
fiow equation for E^, the integration from A = oo to 
A = Aq yields a contribution which does not vanish for 
Aq — > oo, but rather tends to a finite constant. The 
resulting initial condition at A = Aq — > oo reads 



[1/2. 
. 



The right-hand side of Eq. (39) tends to zero fast enough 
so that no such additional contribution is built up for the 

flowing interaction and U^" = U . 

Neglecting the leads, at T = 0, and in the limit of 
a single weak or strong impurity the coupled differen- 
tial equations (39) and (41) can be solved analytically.^^ 
The results are consistent with what is expected from the 
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LSGM. For general parameters we solve the flow equa- 
tions numerically. Wc have developed an algorithm"''® 
which for r > approximately scales with N /T . For fi- 
nite T systems of 10'^ lattice sites are considered, roughly 
corresponding to the length of quantum wires accessible 
to transport experiments. Since the interacting wire con- 
sists of N lattice sites the energy scale 



3JV 



(42) 



forms a lower bound for any temperature scaling of the 
conductance with an interaction dependent exponent (see 
Sect. V) and G(T) quickly saturates for T ^ 6^- There- 
fore only temperatures of the order of the band width 
down to T = 10~^ are relevant and have been stud- 
ied. The Fermi velocity is given by vp = 2sinfcir, with 
the Fermi momentum kp = nn. We note that the sat- 
uration of T,^ for A ^ T or A ^ 5n sets in "auto- 
matically" , in contrast to more intuitive renormalization 
group schemes in which the flow of the considered quanti- 
ties is stopped "by hand" by replacing A — > T or A — > ^jv, 
respectively. ''''^'^'^^'^^ 

At the end of the flow (A = 0), S^^*^ presents an ap- 
proximation for the exact self-energy and will be denoted 
by S in what follows. The fRG is set up in the grand 
canonical ensemble with fixed chemical potential /U. As 
we want to compare our results with those obtained for 
the LSGM considering the canonical ensemble with fixed 
density n we first determine the parameter z^(n, U) such 
that at T = the density on the sites 1 to A'^ acquires the 
desired value n. In a second step the global chemical po- 
tential /i(n, U,T) is determined such that the density of 
the entire system (interacting wire and leads) remains n 
for all T. Due to particle-hole symmetry in the half-filled 
case we have v{l/2, U) = 1/2 and /i(l/2, U,T) = 0. 

Typical results of the A^-dependent T,jj and ^jj+i at 
T = for a single impurity located at site jo and U > 
were presented in Figs. 5 and 6 of Ref. 36. Both matrix 
elements show long-ranged spatial oscillations induced by 
the impurity with an amplitude which asymptotically de- 
cays as 1/1 j — jo|-^^ For finite T this power-law decay of 
the self-energy induced by a single impurity is cut off at 
a scale cx 1/T beyond which the oscillatory part of Sj^ 
and decays exponentially. This is exemplified in 

Fig. 3 in which |AI]j^j+i| = — So| is shown for 

different T and as a function of j — jo- The average value 
of Sjj+i away from the impurity site, which is related 
to the interaction-induced broadening of the band, is de- 
noted by So- In the left panel of Fig. 3, |AEj,j-|-i| is 
shown on a log-log scale where straight lines are power- 
laws, while in the right panel it is presented on a linear- 
log scale where straight lines represent an exponential 
decay. The self-energy was calculated for the case of a 
single hopping impurity with t' = 0.1 located in the mid- 
dle of the interacting wire with A^ = 10^ sites. The other 
parameters are U = 1 and n = 1/2. For a hopping impu- 
rity and at half filling the diagonal part of the self-energy 
vanishes. 



\ ■ 


, 1 1 1 


,".1,1, 



10'' 
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10° 10^ 10^ 10^ 1000 2000 3000 
j-Jo 

FIG. 3: Decay of the oscillatory part of the off-diagonal ma- 
trix element of the self-energy away from a single hopping 
impurity at bond jo, jo + 1- Results for t' = 0.1, jo = 5000, 
N = 10^, U = l,n = l/2 and different temperatures 
T = 10"^ (solid line), T = 10"^ (dotted line), T = 10"^ 
(dashed line), and T = lO"** (dashed-dotted line) are pre- 
sented. The left panel shows the data on a log-log scale, the 
right panel on a linear-log scale. For comparison the left panel 
contains a power- law (j — jo)~^ (thin solid line). 



From E we obtain an approximation for the conduc- 
tance G(T) through our interacting system. We use the 
fact that it is consistent to neglect vertex corrections 
to the Landauer-Biittiker formula when the approximate 
self-energy is frequency independent.^^ This is the case in 
our fRG scheme as discussed in the Appendix. Therefore 
using Eq. (16) does not imply any further approxima- 
tion. To obtain G{T) we have to calculate the effective, 
T- and (^Ar-dcpendont transmission \t(e,T, Sn)\'^ (where 
the dependences on e, T, and Sn are from now on writ- 
ten explicitly) which is expressed in terms of the one- 
particle Green function between sites 1 and A'' [see Eq. 
(15)] and can be calculated by inverting the tridiagonal 
matrix Q — E. This way wc have reduced the many-body 
problem to a single-particle scattering problem off an ef- 
fective impurity potential on the sites ItoN given by the 
T- and (57v-depcndent self-energy. For the detailed under- 
standing of our findings for the single- and double-barrier 
case we use the one-particle scattering theory results de- 
rived in Sect. III. 

A comparison of the fRG results for a single impurity at 
T = to accurate DMRG results and exact scaling prop- 
erties shows that the error due to our various approxima- 
tions is small for interactions such that 1/2 < AT < 1.36,38 
As we discuss next, the same holds at finite T. 



V. RESULTS 
A. Single impurity 

In this section wc report our results for transport 
through a single impurity, and later through a double 
barrier. In the single impurity case we mainly consider 
site impurities; we have obtained similar results for hop- 
ping impurities. We first consider impurities placed suf- 
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ficicntly away from the contact regions and then briefly 
discuss the dependence of the conductance on the impu- 
rity position jo- 

Fig. 4 shows typical fRG results for the T dependence 
of the conductance in the case of a single site impurity of 
strength V. The parameters are U = 0.5, n = 1/2, and 
A'^ = 10''. For T > i? all curves show a 1/T scaling. For 
these temperatures the derivative of the Fermi function in 
Eq. (16) only varies very weakly with e for the relevant 
energies —B/2 < e < B/2 but decreases with increas- 
ing temperature as 1/T. At the same time the integral 
over the band energies of T, ^^v)!^ is only weakly T- 
dependent leading to the observed 1/T behavior. For 
a strong impurity V = 10 and Jjv < T < S, G{T) 
nicely follows a power-law with an exponent (3s{U, n) in- 
dicated in Fig. 4 by the dotted line. More specifically 
we find for all U and n that power-law scaling docs not 
set in before T < B/40. It holds down to T « 5jv/2 
3 • 10~* for N = 10"^) beyond which the conductance 
saturates. The exact value of this lower bound depends 
on the impurity position jo (see below). We here con- 
sider jo ~ N/2. For an intermediate impurity V = 1 
the slope of the data (on a log-log scale) tends towards 
Ps but is still significantly away from it when saturation 
sets in at 6^/2. The slow change of the slope seen in 
Fig. 4 is a general feature of intermediate V. This obser- 
vation is of relevance for the analysis of noisy experimen- 
tal data usually restricted to temperatures within one to 
two orders of magnitude. ''^■''^ Under these conditions a 
transient regime might easily incorrectly be identified as 
the asymptotic power-law regime leading to an exponent 
which is too small. For a weak impurity G(T) stays very 
close to e^/h for all T ^B. 




FIG. 4: Temperature dependence of the conductance for a 
single site impurity with different V as indicated in the legend 
and U = 0.5, n = 1/2, iV = lO". The dotted line shows 
a power-law with an exponent /3s obtained from fitting the 
V = 1Q data. 

In Fig. 5 the U -dependence of the strong- impurity ex- 
ponent ps is shown for n = 1/2 (circles) and n = 1/4 

(squares). We also present results for 17 < 0, where 
G(T) increases for decreasing T and < 0. In this 
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FIG. 5: Dependence of the strong-impurity exponent /3s on 
the interaction U. For comparison the LSGM prediction 2a b 
(solid line) and the leading order (in U) result 2ag (dashed 
line) are shown. 

case G{T) cx T'^^ does not present the asymptotic low- 
energy behavior since for T ^ 0, G(T) is tending to- 
wards e^/h. The larger \U\, the faster G{T) increases. 
At the same time the temperature regime over which a 
power-law exponent can be read off shrinks even if V is 
chosen to be very large. An exponent can be reliably 
read off for U > —0.5 only. From the studies of the 
LSGM^'^° we expect that the strong-impurity exponent 
is2aB = 2(l//^-l). Forn= 1/2, if is given in Eq. (10) 
while away from half filling it can be calculated solving 
the Bethe ansatz integral equations numerically. For 
comparison, 2aB (solid line) and the leading-?/ behavior 

2a^ = 2 m^-^o^i''^^^) (43) 
2TrVF 

(dashed line) are presented in Fig. 5. The exponent Pg 
turns out to be at least correct to leading order in U. Fur- 
thermore, even for intermediate U ^ 1.5, it is still very 
close to 2aB in contrast to the leading-order contribution 
2a^. 

For weak impurities G(r) is predicted to follow a 
power-law scaling 1 - G{T)/{e^/h) cx T'^i^-i) ^j^j^j^ 
holds as long as the correction to perfect conductance 
stays small.^'^° For if - 1 > 0, i.e. ?7 < 0, this scal- 
ing presents the asymptotic low-energy behavior. In 
Fig. 6, 1 - G{T)/{e^/h) is shown for V = 0.01, n = 1/2, 
N = 10"*, and different U. We can read off an exponent 
l3yj{U,n). In Fig. 7, (3^, for n = 1/2 and n = 1/4 is 
compared to 2{K — 1) and the leading-order term — 2q^. 
Note that to leading order 2aB = 2{1 — K), which will 
become important in the comparison of the fRG to the 
"leading-log" method. Our weak-impurity exponent (3^ 
agrees to first order in U with 2(K — 1), but is closer to 
the exact result than — 2a^. 

In the solution of the LSGM the strong- and weak- 
impurity exponents characterizing the scaling of G can 
both be expressed in terms of K. Identifying 13s = 
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is only fulfilled approximately. 



FIG. 6: Weak-impurity behavior of 1 - G{T)/{e^/h) for V = 
0.01, n = 1/2, N = 10* and different U as indicated in the 
legend. 




FIG. 7: Dependence of the weak-impurity exponent /3w on the 
interaction U. For comparison the LSGM prediction 2{K — 1) 
(soUd lino) and the leading order (in U) result —2aB (dashed 
line) are shown. 



2(1/ Ks - 1) and (3^ = 2{K^ - 1) we obtain two £RG 
approximations Ks and for K. In Fig. 8 the rela- 
tive error of the fRG approximations to the exact TLL 
parameter AK/ K = \K — Kg/^ | /K is presented for 
n — 1/2 and n = 1/4. From the figure it is apparent 
that AK (X U^. In our approximation scheme terms of 
order are only partially included^^ and we thus can- 
not expect agreement to higher order. The error AK/K 
depends on n. The half-filled case exhibits the largest 
deviation from the exact result. This can be understood 
from the flow of U^. For fixed C/ > as a function of n, 
(jjA=o _ uyu becomes maximal at n = 1/2, where the 
nearest-neighbor interaction increases by a few 10 per- 
cent (for an example at T = and U = 1 see Fig. 4 of 
Ref. 36). In contrast for smaller fillings, e.g. n = 1/4, 
even decreases during the flow. This explains why our 
weak coupling method works particularly well for small 
fillings. The two approximations Kg and Kyj differ to or- 
der C/^ which shows that within our approach the TLL re- 
lation between the strong- and weak-impurity exponents 
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FIG. 8: Relative error of the fRG approximations for the 
TLL parameter from the strong- l3s and weak- /3w impurity 
exponents at n = 1/2 and n = 1/4. 

In the LSGM with repulsive interaction even a weak 
bare impurity induces a power-law suppression of G with 
exponent 2aB in the limit where all energy scales are 
sent to zero. For weak bare impurities and 1/2 < K < I, 
the crossover scale between weak- and strong-impurity 
behavior goes as F^/C^"^) and becomes very small.^'^'^° 
For the lattice model considered here we have observed 
a similarly small scale. As a consequence, even for the 
fairly large system size of N = 10"* sites and the large 
range of temperatures we can treat, it is impossible to 
directly demonstrate the full crossover for a single set 
of parameters. In the LSGM it can be shown indirectly 
using a one-parameter scaling ansatz 



G = -j^Gk{x) 



x=[T/To{U,n,V)] 



K-l 



(44) 



with a non-universal scale To{U, n, l/).6.8,io pgj. appropri- 
ately chosen Tq the G{T) curves for different V (but fixed 
K) can be collapsed onto the _ft'-dependent scaling func- 
tion Gk{x). It has the limiting behavior Gk{x) (x l — x"^ 
for X — » and Gk{x) oc x~'^I^ for a; — » oo. We ear- 
lier demonstrated one-parameter scaling for the micro- 
scopic lattice model considered within our approxima- 
tion scheme replacing T in the above ansatz Eq. (44) by 
To set up the scaling we have to decide whether 
we take Kg or K^ in the definition of the variable x. For 
a comparison to the analytically known K = 1/2 scal- 
ing function of the LSGM it is advantageous to use Kg, 
since for the corresponding U \U = 2 iov n = 1/2; see 
Eq. (10)] Kg is much closer to the exact K than Kw 
(see Fig. 8). This is what we did in Ref. 38. This choice 
leads to a small deviation from the 1 — x^ behavior at 
small X. For small U the difference in the approximate 
TLL parameters is small and we here use Kyj instead. 
In the limit x — > it implies an exact 1 — x^ behav- 
ior of the fRG approximation for the scaling function. 
A one-parameter scaling plot is shown in Fig. 9. The 
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different symbols stand for different V G [0.01,30]. For 
each V, G{T) was calculated for several T e [10^^. 10"^]. 
The open symbols represent results obtained at half fill- 
ing for U = 0.5 which gives = 0.85. The same K^j 
we find for quarter-filling and U = 0.851. Data for the 
latter parameters arc shown as filled symbols. The col- 
lapse of the results for n = 1/2 and n = 1/4 exemplifies 
that the scaling function depends on U and n only via 
the TLL parameter K as predicted for the LSGM. It is 
very remarkable that our exponent 2(5 g/ (5^ of the large 
a;-behavior agrees to leading order with the correspond- 
ing exponent of the LSGM. This is a significantly 
stronger result than the observation that fis and 2aB as 
well as fiw and 2{K — 1) agree to order U. Considering 
site (hopping) impurities which are extended over more 
than one site (bond) as well as impurities which are given 
by a combination of modified on-site energies and hop- 
ping matrix elements we have verified that Gk is also 
independent of the details of the local impurity. 




10 10 10 

x=[T/TJU,n,V)]^«/^ 



FIG. 9: One-parameter scaling plot of the conductance. Open 
symbols represent results obtained for U = 0.5, n = 1/2, and 
different T and V , while filled symbols were calculated for 
U = 0.851, n = 1/4. Both pairs of U and n lead to the same 
Ky, = 0.85. The solid lines indicate the asymptotic behavior 
for small and large x. 

A comparison of the results presented here for fixed 
5n and varying T and the results of Rcf. 38 obtained for 
T = and different 5n shows that in the single-impurity 
case and for smooth contacts 5jv and T present equiv- 
alent scaling variables. In alternative RG schemes this 
equivalence is instead assumed. Note that the equiva- 
lence no longer holds if either non-perfect contacts^^'^^ 
or resonant tunneling is considered (for the latter, see 
below). 

In Ref. 9 an alternative fermionic RG method was in- 
troduced to study an effective low-energy model with a 
single barrier in the limit of weak interaction. In this ap- 
proach an RG equation is set up to resum the "leading- 
log" divergences of the effective transmission at the chem- 
ical potential. The results obtained using the "leading- 
log" method only partly agree to ours. While in our 
approach the one-parameter scaling function depends on 



K and the exponent —2/K of the large .x-behavior is 
reproduced correctly to order [/, the one-parameter scal- 
ing function of the "leading-log" approach turns out to 
be the non- interacting {K =1) function independently of 
the interaction strength chosen.^ This method thus does 
not capture all leading-order interaction effects. 

This shortcoming of the "leading-log" method can be 
traced back to the analytical form of the effective trans- 
mission at scale A 



|tu|-|A/Ao|^"^ 
|roP + |ioP|A/Ao|2«^ 



(45) 



as it follows from the RG equations.® The impurity pa- 
rameters (height and width) enter via the non-interacting 
transmission and refiection probabilities at the chemical 
potential |to|^ and |ro|^, and Aq denotes a non- universal 
energy scale. As above the exponent a^g stands for the 
leading-order term oi as = 1/K — 1. This equation is 
derived under the assumption that all energy scales, i.e. 
T, 5n, and the energy of the incoming particle e (mea- 
sured relative to the chemical potential) are set to zero. 
The dependence on T, Sn, or s is then determined by 
replacing A by one of these variables. For T = 0, at the 
chemical potential, but for a finite system this gives 



\t{SN)\' 



\to\HSN/So 



,2a' 



\ro\^ + \to\HSN/So) 



2a' 



(46) 



with the non-universal scale Sq. Applying Eq. (16), which 
for T = reduces to G = e^|f(/i)p//i, and introducing 

X = Iro/foKf^Ar/^o)^"^; where — a'^ is the leading-order 
term oi K — 1, the scaling function of the conductance 
turns out to be the non-interacting one, Gk=i = 1/(1 + 
■T^). Eq. (46) can be compared to our Eq. (24) which 
applies to the situation of a single site impurity {Vg = V). 
Because of the above-mentioned equivalence of the T and 
5n scaling it is sufficient to consider T = and the length 
of the interacting wire as the variable. Without loss of 
generality wc here focus on n = 1/2 with = 0. For 
ii = U, T — 0, and at the chemical potential, i.e. for 
£ = 0, Eq. (24) reads 



|i(0,0,<5,v)|' 



4A2(0,0,(5^v) 



y 4- 217(0, 0, (5jv) -I- 4A2 (0, 0, (5 



}N) 



Because of the left-right symmetry we here and in the 
following suppress the index l/r. During the fiow the 
potential V on the impurity site gets renormalized. At 
A = it is denoted by V . Eq. (47) holds for a more 
general impurity in which also the hopping amplitude 
ti = tr [see the discussion following Eq. (9)] to the left and 
right of the lattice site with onsite energy V is different 
from one.^® A and fl can be obtained numerically from 
the self-energy at A = as explained in Sect. III. For 
F — > 0, which implies V ^ because of the particle- 
hole-symmetry for V = 0, and sufficiently large N we 
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find 

V + 20(0, 0, ^jv) oc ^^"/^ , A(0, 0, ^jv) = const. (48) 
while for y » B 

F + 20(0, 0, (Jjv) = const. , A(0, 0, (5jv) (5^"^^ . (49) 

Instead of the one exponent a'g of the "leading-log" ap- 
proach we obtain two different exponents characteriz- 
ing the weak- and strong- impurity limits. This leads to 
an interaction-dependent scaling function with a large 
X scaling exponent which is correct to leading order in 
U. In contrast to Eq. (46) the analytical form Eq. (47) 
allows for two different sources of an A^-dcpcndcnce of 
the effective transmission and the conductance. We thus 
conclude that reducing the number of flowing coupling 
constants to a single one — the effective transmission at 
the chemical potential — leads to an analytic form of 
the conductance which is too restricted to cover all im- 
portant leading-order interaction effects. This shows that 
our fRG approach goes beyond the "leading-log" method. 



1.5 



<l 




|V+2n|/2 

FIG. 10: Flow diagram of the conductance in the plane de- 
fined by the real and imaginary part of the auxiliary Green 
function. The direction of the flow for Sn ^ and ?7 > is 
indicated by the arrows. For U < it is reversed. The axis 
are lines of fixed points. A-a^xis: line of "perfect chain" fixed 
points with G = /h. {V + 2f2)-axis: line of "open chain" 
fixed points with G = 0. 

Eq. (47) can be used to further visualize the scaling 

of the conductance. Within OTir approximation the con- 
ductance can be expressed in terms of the two scale- and 
parameter-dependent functions {V + 20)/2 and A > 0. 
Fig. 10 shows a parameter plot in the ([V" + 20]/2, A) 
plane obtained at e = T = and for varying N . Each 
curve was calculated for a fixed V and ti = tr &t U = 1 
and half filling. A similar plot can be obtained for a large 
fixed iV calculating {V + 20) /2 and A for different A > 
considering T,^ as an effective scattering potential. With 
increasing N, i.e. decreasing energy scale, the "flow" ap- 
proximately follows a section of a circle centered around 
the origin, with a radius which depends on the initial im- 
purity parameters. For ?7 > the "flow" is clock-wise. 



For attractive interaction the direction is reversed. The 
axis with {V + 20)/2 = can be considered as a line 
of "perfect-chain" fixed points with |tp = 1, while the 
axis with A = represents the line of "open-chain" fixed 
points with \t\'^ = 0. In accordance with Eqs. (48) and 
(49), close to the two lines of fixed points the "flow" is 
perpendicular to the axis. 




FIG. 11: The conductance G(T) for a site impurity with V = 
10, U = I, n = 1/2, AT = 10" -I- 1, and different impurity 
positions jo. 

The above results are generic as long as the impu- 
rity is placed sufficiently away from the contact regions. 
A richer behavior is found if the impurity is positioned 
closer to the smooth contacts. Because of symmetry we 
only have to consider the case in which the impurity is 
moved towards the left contact. Then the scale 

^.0 = ^ (50) 

becomes important. -'^^ It sets the lower bound for the 
power-law scaling with the exponents discussed above. 
For T « i5j|j we find a crossover to a power-law scaling 
with different exponents. This is exemplified in Fig. 11 
for a strong impurity with V = 10, A'' = 10'* + 1,U =1 
and three impurity positions jo = 5001 (solid line). 156 
(dashed line), 56 (dashed-dotted line). For jo = 156 at 
T w the scaling with exponent (is crosses over to a 
scaling with exponent Psl'^- This coincides with the pre- 
diction obtained within a field theoretical model where 
the exponent 2a b crosses over to half its value as-^^ For 
impurity positions in the contact region, i.e. jo = 56, we 
only observe the exponent /3s /2 that describes the tun- 
neling between the non-interacting leads and the TLL. 
The curve for jo = 56 does not follow the other two 
curves even at higher temperatures as compared to the 
bulk value the interaction at jo is reduced by a factor of 
two. 

Similar behavior is found for the scaling of 1 — 

G{T) / [e^ /h) in the limit of a weak impurity. In this 
case the exponent predicted for Jjv < T < is 
2{K - 1) / (K + I) and \s thus different from half the 
exponent 2{K - 1) found for 5j^^ < T < B.^^ The ex- 
ponent we find agrees to leading order in U with this 
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prediction. We note that even though (3^/2 agrees with 
2{K — 1)/{K+1) to order U our exponent in this regime 
is different from /3u,/2. 



B. Resonant tunneling 

In this section we only investigate repulsive interac- 
tions U > Q and focus on n = 1/2, with /x = 0, as the 
results for other fillings do not differ qualitatively. Wc 
furthermore only consider situation in which the quan- 
tum dot is placed sufficiently away from the contact re- 
gions. 



infinitely sharp. It is important to note that for all pa- 
rameter sets tested we found agreement between twice 
the scaling exponent of ct at T = and the weak-single- 
impurity exponent up to our numerical accuracy. The 
difference of cr for the different resonance peaks at a fixed 
N (see Fig. 12) is a band effect. Off resonance, G vanishes 
as 



(52) 



1. T — 

We start our investigation of resonant tunneling con- 
sidering T = 0. Fig. 12 shows typical results for the gate- 
voltage Vg dependence of the conductance for symmetric 
barriers. The quantum dot contains Nd = 6 sites and 
is separated from the rest of the interacting wire by two 
site impurities of intermediate height Vi/r = 2. In the (51) to leading order in U . 
figure the A^-indepcndent [/ = result (dashed-dotted 
line) is compared to the conductance of fRG calculations 
with U = 1 for two different A''. For not-too-weak barri- 
ers the number of perfect resonances (with peak conduc- 
tance e^/h) at positions Vg coincides with the number 
Nd of lattice sites which form the dot. For intermediate 
barriers and [/ = the peaks still overlap. For increasing 
Vijr the width of each peak decreases until well-separated 
resonances are formed. Also for {/ > 0, G shows perfect 
resonances with peak conductance Gp = jh. This be- 
havior can be understood from Eq. (27), which for ?7 > 
and sufficiently large N can be used even for weak to in- 
termediate barriers as the resonance peaks do no longer 
overlap. This was the necessary condition for the deriva- 
tion of Eq. (27). At VI we find + 217(")(0, 0, 8n) = 0. 

Regardless of the A^-dependence of A(")(0, 0, Jjv) this 
leads to \t{0,0,6N)\'^ = 1. At T = only e = con- 
tributes in the Landauer-Biittiker formula (16), leading 
to G = e'^/h. Remarkably, A(") (0, 0, ^at) at does not 

/2 

follow the scaling Sj^ of the spectral weight close to a 
single impurity^ but instead saturates at small Sn- The 
imaginary part of the exact Green function, which can 
be calculated from the auxiliary Green function using 
Eq. (28), behaves similarly. Increasing U at fixed N re- 
duces the overlap between the resonances by increasing 
the energy difference between the peaks but even more 
importantly by reducing the width of each peak. The 
half- width a depends on N. We find that the width of 
all peaks scales to zero as 



with the strong-single-impurity exponent (3s (again up to 

the numerical accuracy). Using Eq. (27) this behavior 

can be traced back to a A("^(0,0,(5]v) cx: S^^^ scaling 
while Sa + 2f2(")(0, 0, Jat) goes to a constant different 
from 0. Off resonance the double barrier acts as a single 
impurity of strength AVg = Vg — Vg With this insight 
and using the one-parameter scaling ansatz Eq. (44) the 
A^-depcndence of cr, Eq. (51) can be explained as follows: 
For small \AVg\ the non-universal scale (here i5o) goes as 

1/2 this leads to xi/2 = 

thus find cr 



With G^(a;i/2) 
AVg^^^^ I ) with a AV^-independent constant c. We 
2|AV"J^/^^| oc which agrees with Eq. 




cr OC ^ 



N 



(51) 



with the wcak-single-impurity exponent (3^ as discussed 
in the last section. For N ^ oo the resonances are thus 



FIG. 12: Gatc-voltage dependence of the conductance for a 
symmetric double barrier at T = 0. The parameters are No = 
6 and Vi/r = 2. 

For asymmetric barriers and U = the conductance 
still displays resonances but with peak values Gp which 
are smaller than e'^ /h. The same holds for J7 > as 
is shown in Fig. 13 for hopping barriers with ti = 0.1, 
tr = 0.3, U = 0.5, and Nd = 10. We set the in- 
teraction across the barrier bonds to 0. The number 
of resonance peaks is Nd- Here we only show a single 
peak. For hopping impurities as barriers, G is symmet- 
ric around Vg = (for n = 1/2). With increasing U 
and for fixed N, Gp decreases. The same holds for fixed 
U and increasing N. For sufficiently large N we find 
that Eq. (52) holds for the off- and on-resonance scal- 
ing, i.e. on small scales also at VJ an asymmetric double 
barrier acts as a single impurity. As a consequence for 
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FIG. 13: Gate voltage dependence of the conductance for an 
asymmetric double barrier at T = 0. The parameters arc 
U = 0.5, Nd = 10, ti = 0.1, and tr = 0.3. Only gate voltages 
close to Vg are shown. Note the linear-log scale. 



N ^ oo the half-width cr of the peaks approaches a con- 
stant. Again Eq. (27) can be used to further analyze 
G. At V^, Ea + (0, 0, 5n) + f^r"' (0, 0,6n) = 0. For 
tr < ti < 1 and 5n — > 0, Ar{0,0, 6n) vanishes as 6^^^^ 



-A/2 



N 



These results corn- 



while A;(0, 0, (5 at) increases as S 

bine to Eq. (52). Applying Eq. (28) one can show that 
in contrast to the asymmetric behavior of the imaginary 

part Aj/j, of the auxiliary Green fimction, the spectral 
function of the exact Green function on the sites to the 
right and the left of the barriers both scale to zero as 

3 /2 

S]^ . The off-resonance behavior can be explained sim- 
ilarly as for symmetric barriers: Ai/j.{0,0, 6n) oc S^^"^ 

while Ea + r^|"^(0, 0, (Jat) -I- fli°'\o,0, S^) goes to a con- 
stant different from 0. We obtain all the above results 
for both hopping and site impurities as barriers. 

Within an extension of the LSGM to the two-barrier 
case resonant tunneling on asymptotically small scales 
has been discussed in Ref. 6. The results of that work 
agree with ours. As for the single-impurity case we re- 
produce the exponents of the power-laws (51) and (52) 
to leading order in U, i.e. the exponents of the extended 
LSGM are given hj 1 — K and 2aB, respectively. In 
the low-energy limit, the conductance exhibits the same 
scaling behavior as a function of T or Sn, but for higher 
temperatures the T dependence of G at fixed A'' is much 
richer than the Sn scaling, showing non-monotonic be- 
havior and power-laws with different exponents in a va- 
riety of temperature regimes. 



2. T > 

Shortly after the discussion of resonant tunneling on 
asymptotically small scales the conductance for interme- 
diate to high energies was investigated within a field- 
theoretical model using different approximate methods. 



For symmetric barriers and applying second-order per- 
turbation theory in the barrier height the T > devia- 
tions of Gp{T) (at resonance) from e^/h were determined 
to scale as T"^^.^^ At further increasing T a regime of 
uncorrelated sequential tunneling (UST) was predicted 
based on a perturbative analysis in the inverse barrier 
height, where Gp{T) oc T'^b-i and a{T) oc T.^^'^o j,^^^ 
regime is bounded from above by the (non-interacting) 
level spacing of the dot Sn^ = ttvf/Njj. Within the 
same perturbative approach and for d^j^ <C T <C -B, 
Gp{T) increases as T^"^ for increasing T}"^ In contradic- 
tion to the temperature dependence following from UST, 
in recent transport experiments on carbon nanotubes a 
suppression of Gp with decreasing T was observed for 
T ^ ^jvd-^^ Using another approximation scheme, "cor- 
related sequential tunneling" (CST) was predicted to re- 
place UST.21.27 It leads to Gp cx t^»b-\ Although a 
simplified field-theoretical model of spinless fermions was 
used and contacts as well as leads were ignored the CST 
scaling was argued to be consistent with the experimental 
data. This consideration has stimulated a number of the- 
oretical works. 2^ In particular, Quantum Monte Carlo 
(QMC) results^^ for Gp{T) were interpreted to be consis- 
tent with CST for T ^ 6nj^ . Resonant tunneling was also 
investigated using the "leading-log" method. ^^'^^ Within 
this approach no indications of a CST regime were found. 
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FIG. 14: Schematic plot of the different regimes for the scaling 
of the conductance Gp at resonance for symmetric barriers 
found using our approximation scheme. The upper left panel 
is for dots with a large number of lattice sites Nd and high 
barriers, the upper right for small Nd and high barriers, the 
lower left for large No and low barriers, and the lower right 
for small Nd and low barriers. The relevant energy scales are 
the band width B, the U = level spacing of the dot 5nd, 
the lower bound of the UST regime T* [see Eq. (53)], and the 
infrared cutoff <5jv associated with the number of lattice sites 
forming the interacting wire. 

Being applicable for any barrier height, but only for 
interactions such that 1/2 < if < 1, our fRG based 
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method is complementary to the approaches of Refs. 17 
and 20. It has the advantage that all relevant tempera- 
ture regimes can be investigated within one approxima- 
tion scheme (see Ref. 40 and below). 

In Fig. 14 our findings for Gp{T) and a fixed N pre- 
sented in Ref. 40 are summarized schematically in the 
four limiting cases of small and large quantum dots sep- 
arated from the rest of the interacting wire by high and 
low barriers. The behavior found for scales of the 
order of B and larger is the same band effect as discussed 
above for the single- impurity case. 

For large dots and high barriers (upper left panel of 
Fig. 14) the T^^ behavior is followed by a T^^" scaling, 
with the strong-single-impurity exponent (3s, down to 
temperatures slightly less than 6nj^. According to Eq. 
(32) in this parameter and temperature regime Gp{T) is 
obtained by adding the resistances of the two individual 
barriers (Kirchhoff's law). Both follow the single impu- 
rity scaling T^' which explains that the same exponent 
is found in resonant tunneling. For even lower T, Gp fol- 
lows a power-law with (to leading order in U) the UST 
exponent /3s /2 — 1 down to a scale T*. The crossover 
from the regime in which Kirchhoff's law can be applied 
to the UST regime is relatively sharp and takes typi- 
cally half an order of magnitude (see Fig. 2 of Ref. 40). 
For T <C 5jvd the width of —df/de is smaller than 
and only a single resonance peak (the one around 0) of 
|t(£, T, (5jv)p contributes to the integral in Eq. (16). The 
width A^") of this peak vanishes as t^T^'^'^/Nd leading 
to Gp{T) oc T^»/2-i as discussed in Sect. III. In the case 
of hopping barriers r stands for ti/j. while it is propor- 
tional to l/Vi/r for site barriers. The lower bound of this 
scaling regime is reached when T is equal to A^"', i.e. at 

T* oc {t^/NdY^^^-^^^^^ . (53) 

Using Eq. (16) it is easy to see that for T* <C T ^ Sno 
the half-width a of the resonance in G{T) as a function 
of Vg scales like T. These results hold for all cases in 
which we find a UST regime, i.e. for the other dot pa- 
rameters discussed next, but also for asymmetric barriers 
and the off-resonance scaling investigated further below. 
For large dots, high barriers, and the values of N used 
here, T* is very close to the scale S^. Below this scale 
no interaction-dependent "universal" behavior can be ob- 
served. 

For small dots and high barriers (upper right panel of 
Fig. 14) Snd s-iid B are very close to each other and no 
regime with T'-^" scaling is developed. Due to the factor 
1/Nd in Eq. (53), for small dots T* is larger than for 
large A''^. For small dots a UST regime with exponent 
/3s /2 — 1 is not as nicely developed as for large Njj (see 
Fig. 2 of Ref. 40 and Figs. 16, 17, and 18 below). 

For large dots but low barriers (lower left panel of Fig. 
14) T* moves to higher energies because of the factor 
T in Eq. (53). Between B and 6no^ Gp decreases with 
decreasing temperature. As is discussed in Ref. 40 this 
decrease is not described by a power-law. We thus denote 
it as non-universal behavior. 
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FIG. 15: Effective exponent (i.e. logarithmic derivative) of 
1 - Gp{T)/{e^/h) for a small dot with No = 1, U = 1, N = 
W^, and intermediate to weak hopping barriers ti/r- Dashed- 
dotted line: 2Ku,. 

For small Njy and low barriers (lower right panel of Fig. 
14) 6n and T* are sufficiently separated so that we can 

observe an additional temperature regime with power- 
law scaling of Gp and an interaction-dependent exponent 
as exemplified in Fig. 15 for [/ = 1, A = 10^, = 1, 
and the barrier heights ti/r = 0.5 as well as t;/^ = 0.7. 
We find that for Jat < T < T*, Gp{T) approaches e^/h 
as 

1 - Gp{T)/{e^/h) oc T2^» , (54) 

with the fRG approximation Jfu, for K obtained from 

f3,u] (sec Sect. VA). In this temperature regime and 
for energies e which contribute to the integral in Eq. 
(16), |t(e,T,5jv)P does not depend on T. We find that 
for (^Tv < e < T*, e - e„ - 2n(") (e, 0, (^w) oc e and 
A^"\s,0,Sjm) cx e^'^"'/^. Evaluating the e integral for 
(^w < T < T* leads to Eq. (54). For low barriers, smaU 
Nd, and T* < T < ^jv^ we no longer find the UST 
power-law but instead non-universal behavior. 

For the cases in which a comparison is possible the 
above exponents agree (to leading order in U) with the re- 
sults obtained for a field-theoretical model using pertur- 
bation theory in the barrier height and the inverse barrier 
height. ^^'^^ In addition applying the fRG to our micro- 
scopic model we can describe the complete crossover be- 
tween the temperature regimes with "universal" power- 
laws and cover the parameter regimes with non-universal 
behavior. 

We do not find any indications of the CST exponent 
2aB — 1. Besides the single- and double-barrier problems 
some of us also studied transport of correlated electrons 
through junctions of M Id wires, e.g. Y-junctions with 
M = 3 and X-junctions with M = 4.^" Also in these ge- 
ometries one finds power-law scaling of the conductance. 
For a certain case (Y-j unction with magnetic fiux) the in- 
dependent exponent is known analytically,^^ which we 
again reproduce to leading order in We thus con- 
clude that in all cases of inhomogeneous TLLs we have 
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studied and in which exact resuhs for exponents (from ef- 
fective field-theoretical models) of power-law scaling are 
known we reproduce these exponents at least to leading 
order. We find it thus very likely that the occurrence of 
the CST regime with exponent 2aB — 1 is an artifact of 
the approximations used in Refs. 21 and 27. 

As Fig. 3 of Ref. 40 shows, the non-universal behavior 
discussed above might quite easily be identified incor- 
rectly as power-law scaling with an exponent, which is 
significantly smaller than j3s ~ 'i^as, i-e. in the vicinity of 
2aB — 1 . This is of particular importance for the interpre- 
tation of noisy QMC data as well as noisy experimental 
results, both typically restricted to a temperature regime 
of one to two orders of magnitude. ^^■''^ In Fig. 4 of Ref. 
26 a single set of QMC data with if = 0.6, for an inter- 
mediate dot size, and intermediate barriers is presented 
roughly falling into the parameter regime shown in the 
lower left panel of Fig. 14. In Ref. 26 a certain part of the 
temperature regime with decreasing Gp{T) (for decreas- 
ing T) was fitted by a power-law with exponent 2aB — 1 
which was then interpreted to support the occurrence of 
CST. 




FIG. 17: Upper panel: G{T) for symmetric barriers with 
ti/r = 0.2 for ?7 = 1, iV = 10", ATo = 4 at resonance |AFg| = 
(circle), close to a resonance with |AVg| = 0.005 (squares), 
and further away from the resonance with |AVg| = 0.1 (dia- 
monds). Lower panel: Logarithmic derivative of G{T). Solid 
line: fis', dashed line: /3s/2 — 1. 
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FIG. 16: Upper panel: G{T) for symmetric barriers with 
Vi/r = 10 for U = 0.5, N = 10*, No = 100 at reso- 
nance lAVgl = (circle), very close to a resonance with 
\AVg\ = 0.001 (squares), and in a conductance minimum with 



|AV"g| = 0.04 (diamonds), 
tive of G(T). Solid line: /£ 



Lower panel: Logarithmic derivar- 
dashed line: /3s/2 — 1. 



We next discuss the T dependence of the conductance 
away from Vg still considering symmetric barriers. In 
the upper panel of Fig. 16 G{T) is shown for U = 0.5, 

= 10^, a large dot with Nu = 100, high barriers with 
Vi = Vr = 10, and three different |AV^|. Results for the 
resonance closest to Vg = are shown. We have verified 
that similar behavior can be found for the other reso- 
nances. The circles show the conductance at |AV^| = 0. 
The peak conductance Gp{T) displays the power-laws 
discussed above. This is best seen in the lower panel 
in which the logarithmic derivative of the conductance is 
shown. Close to the resonance (squares) for decreasing 
T, G{T) follows Gp(T) into the UST regime before the 



conductance at |AV^| ^ decreases while Gp(T) tends 
towards e^/h. For |AVg| ^ the peak in \t{e,T,SN)\'^ is 
centered around AVg. In the UST regime its width still 
scales as A^") cx t'^T'^"/'^ /No (sec Eq. (27) and above). 
As long as |AV^| <C 5nd the curves separate when A^") 
is equal to |AV^|. This leads to the "off-resonance" tem- 
perature scale 



n,<x{ND\AVg\/T^ 



a IPs 



(55) 



at which the off-resonance conductance separates from 
Gp{T). Sufficiently away from the resonance, e.g. in 
the conductance minimum (diamonds), the temperature 
regime in which Kirchhoff 's law can be applied is followed 
by a regime in which the double barrier acts as a strong 
single impurity. This crossover typically takes one order 
of magnitude in T and consists of a drop of G by at least 
two orders of magnitude. In both these regimes G{T) 
follows a power-law with exponent which arises for 
completely different reasons. The origin of the power-law 
for (57Vj3 < T <C -B has been explained above. For these 
temperatures many peaks in \t{e,T,6N)\'^ contribute to 
the e integral in Eq. (16). The single-impurity scaling 
on asymptotically small scales can be understood as fol- 
lows: Off resonance and close to the chemical potential 

\t{e,T,5N)? behaves like [e"^ + T'^f''''^ (for e,T > ^jv)- 
In the evaluation of the integral in Eq. (16) the derivative 
of the Fermi function can be considered as being a con- 
stant (as a function of e) which scales like l/T. The 

remaining integral with integrand (e^ + T"^^^"^"^ runs 
over an energy range of order 2T which then leads to 
G{T) oc Tl^'. For T < 6^/2 a deviation from the de- 
scribed power-law scaling sets in. As an additional ex- 
ample for the off-resonance behavior, in Fig. 17 we show 
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G(T) for U = I, Nd = N = 10^, weak hopping bar- 



riers with t 



l/r 



0.2 (again the interaction across the 



barriers is set to 0) and three different |AV^| 



[ jBnQBmQQnHQeBElfii 
10"^<>O«oO*' 




FIG. 18: Upper panel: Gp{T) for asymmetric barriers with 
U = 0.5, N = 10'^, Nd = 2, and Vi = 7, K = 12.3 (squares) 
as well as Vi = 4, Vr = 13.6 (diamonds). For comparison 
also the symmetric case with Vi = Vr = 10 is shown (circles). 
Lower panel: Logarithmic derivative of Gp(T). Solid line: f3s; 
dashed line: /3s /2 — 1. 
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FIG. 19: Upper panel: Gp{T) for asymmetric barriers with 
Vi = 5, K = 50, U = 0.5, N = 10^', No = 2 (circles) and 
Nd = 20 (squares). Lower panel: Logarithmic derivative of 
Gp{T). Solid line: A; dashed line: /3s/2 - 1. 

We finally describe the T-dependence of the conduc- 
tance for asymmetric barriers at resonance. In Fig. 18 
results are presented for U — 0.5, = 2, N = 10^, 
and different site barriers with weak to intermediate 
asymmetry Vi = 7, Vr = 12.3 (squares) and V; = 4, 
Vr = 13.6 (diamonds). The combination of VJ/^ was 
chosen such that at high temperatures the curves ap- 
proximately collapse. For comparison also the symmet- 
ric case with V; = = 10 is shown (circles). Over a 
wide range of temperatures Gp{T) for asymmetric bar- 
riers shows the same behavior as in the symmetric case. 
The scale Tag on which Gp{T) for Vi Vr starts to devi- 
ate from the behavior of the symmetric case discussed 



above decreases with vanishing asymmetry parameter 
7 = |1^2 _ v;2|/(T/2 _^ addition T^^ decreases 

for increasing No as can be seen in Fig. 19 which shows 
Gp{T) for U = 0.5, = 5, = 50, iV = 10^ and two 
different Nu. Even for a fairly strong asymmetry with 
7 « 0.98 {Vi = 5, Vr = 50) for TV = 10'* we observe 
the asymptotic low-energy behavior Gp{T) oc T^' only 
for very small Nd (see Fig. 19). Already for Nu = 20, 
Tks < Sn- 

The temperature dependence of the conductance for 
asymmetric barriers and Vg ^ Vg can be deduced from 
Gp{T) obtained in the last paragraph (asymmetric bar- 
riers) and G{T) for symmetric barriers at \AVg\ > dis- 
cussed further above. 



VI. SUMMARY AND PERSPECTIVES 

In this paper we presented a method to calculate the 
linear conductance through inhomogeneous TLLs con- 
nected to semi-infinite non-interacting leads. The sys- 
tem is described by a microscopic lattice model of spin- 
less fermions. As applications we studied transport in the 
presence of a single impurity as well as the double-barrier 
problem, allowing for resonant tunneling. The contacts 
are modeled to be smooth, i.e. in the absence of impu- 
rities and for T -C B the conductance is equal to e^/h. 
This requires that the interaction is turned on smoothly 
close to the two contacts. The transport problem was 
previously investigated using effective field-theoretical 
models. In the cases where exact results for such models 
are known we obtained quantitative agreement for weak 
to intermediate interactions with 1/2 < if < 1. At half- 
filling this corresponds to nearest-neighbor interactions 
< U < 2, which covers the whole parameter regime in 
which the model is a TLL (with K < 1). Away from 
half-filling even larger U can be considered. 

Our method captures the expected power-law scaling 
of the conductance for a single weak and strong impurity. 
In addition, we can describe the complete crossover gov- 
erned by a K-dependent one-parameter scaling function. 
For a single impurity we briefly studied attractive inter- 
actions. Also in this case our method provides reliable 
results for not-too-strong interactions with 1 < K < 3/2, 
i.e. at half-filling for -1 < ?7 < 0. 

For resonant tunneling depending on the parameters 
of the quantum dot we find several temperature regimes 
with power-law scaling as well as non-universal behav- 
ior. All these temperature regimes are obtained within 
the same approximation scheme. The crossover between 
the regimes can be studied in detail. For parameters 
for which a comparison is possible our results agree with 
the ones obtained in lowest order perturbation theory in 
the barrier height and inverse barrier height. We do not 
find any indications of a "correlated sequential tunnel- 
ing" regime with the exponent 2aB — 1 predicted from an 
approximate Master-equation approach. 2*'^'' If it would 
be present our method should be able to reveal such a 
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regime with a scaling exponent which is of leading order 
in the interaction. 

Our method is a very flexible tool to investigate models 
of inhomogeneous TLLs. It allows for extensions which 
are required to obtain a more realistic description of ex- 
perimental systems: (i) The spin degree of freedom can 
be included, (ii) Models with more realistic contacts (ex- 
tended bulk and end contacts) can be treated. After 
projecting out the non-interacting leads the microscopic 
details of the contacts enter the fRG flow equations as the 
initial condition of the self-energy. Preliminary results 
show that imperfect leads have a strong effect on G(T) 
even in cases in which the imperfection is fairly weak, 
(iii) More realistic models for the leads can be used. All 
that enters the fRG approach is the local Green function 
of the lead at the contact after disconnecting the lead and 
interacting wire. In addition the transport through junc- 
tions of interacting wires, such as Y- and X-junctions, 
as well as transport through Aharonov-Bohm geometries 
can be described using the fRG. 
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Appendix 

The Kubo formula for the conductance has a contribu- 
tion with bare current operators and one with current- 



vertex corrections. In this appendix we explain why the 
latter vanish in our approximation, consistent with the 
Ward identity. 

We have approximated the full two-particle interaction 
vertex F by a renormalized nearest-neighbor interaction 
which is independent of frequency. Then the flow of the 
right current vertex A/f(a;-|-iO) is of order ui and vanishes 
in the limit of the dc conductance. The approximation 
that F is frequency independent has another consequence: 
by the flow equation, it follows that ImS = 0, that is we 
do not capture inelastic processes at second order in the 
interaction. They could be included in the flow equation 
by retaining the frequency dependence of F. 

However, the fact that the vertex corrections and 
ImS vanish simultaneously shows that our approxima- 
tion is at least consistent with the Ward identity asso- 
ciated with particle number (charge) conservation. The 
global continuity equation for the interacting region is 
dnc/dt + Jr — Jl = 0, where nc is the total particle 
number in the interacting region, and Jl.r are the cur- 
rents entering the system from the left and leaving it 
to the right. This continuity equation implies a Ward 
identity.^'' After analytical continuation to frequencies 
slightly above and below the real axis (the branch cut 
of the current operators and the self-energy), the density 
response drops out and we obtain 



PR(e + iO,e- iO) -PL{e + iO, e - iO) = 2 ImS(e -h iO) , 



where Pl,r are the corrections of the left and right cur- 
rent vertices, using the notation of Ref. 57. Thus, the 
current-vertex corrections are related to the discontinu- 
ity across the real axis of the self-energy. At T = 0, both 
sides vanish exactly because there is no inelastic scat- 
tering, while at finite temperatures, neglecting inelastic 
processes is a consistent approximation. 
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